Dynamical features of Shannon information entropy of bosonic cloud in a tight trap 
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We calculate Shannon information entropy of trapped interacting bosons in both the position 
and momentum spaces, Sr and 5**; respectively. The total entropy maintains the fuctional form 
5* — a + blriN for repulsive bosons. At the noninteracting limit the lower bound of entropic 
uncertainty relation is also satisfied whereas the diverging behavior of Sr and Sk at the critical 
point of collapse for attractive condensate accurately calculates the stability factor. Next we study 
the dynamics of Shannon information entropy with varying interparticle potential. We numerically 
solve the time-dependent Gross-Pitaevskii equation and study the influence of increasing nonlinearity 
in the dynamics of entropy uncertainty relation (EUR). We observe that for small nonlinearity the 
dynamics is regular. With increase in nonlinearity although Shannon entropy shows large variation 
in amplitude of the oscillation, the EUR is maintained throughout time for all cases and it confirms 
its generality. We also study the dynamics in a very tight trap when the condensate becomes highly 
correlated and strongly inhomogeneous. Time evolution of total entropy exhibits aperiodic and 
fluctuating nature in very tight trap. We also calculate Landsberg's order parameter for various 
interaction strengths which supports earlier observation that entropy and order are decoupled. 

PACS numbers: 89.70.Cf, 03.75.Kk 



I. INTRODUCTION 

Statistical correlation and Shannon information en- 
tropics in both the position and momentum spaces are 
the key quantities in the understanding of quantum me- 
chanical systems which measure information regarding 
localisation of the position and momentum distribution. 
An important step in this direction is the entropic uncer- 
tainty relation (EUR). Using the position and momen- 
tum space entropies (Sr and Sk respectively), Bialynicki- 
Birula and Myciclski (BBM) derived a stronger version 
of the Heisenbcrg uncertainty relation [ij. For a three- 
dimensional system the total entropic sum has the form 



S for the total density distribution in different quantum 
systems is proposed which takes the form Q 



S^a + b InN: 



(2) 



S = Sr + Sk>3{l+ In (tt)) ~ 6.434. 



(1) 



This means that the conjugate position and momentum 
space information entropies Sr and Sk maintain an in- 
verse relationship with each other. It signifies that when 
a system is strongly localised in position space, the corre- 
sponding information entropy and uncertainty in position 
space decrease. The corresponding momentum distribu- 
tion becomes delocalised, i.e. the information entropy 
and uncertainty in momentum space increase. It is im- 
portant to note that BBM inequality (Eq.(l)) presents a 
lower bound and the equality is maintained for a gaussian 
wave function. 

In the past, work has been devoted to the study of the 
Shannon entropies in various quantum mechanical sys- 
tems like nuclei, atomic clusters, fermionic and bosonic 
systems 0-Ql- The universal property of total entropy 



a and b are the parameters which depend on the given 
systems. Information theory in a system of correlated 
bosons in a trap and various features of information en- 
tropy of confined atoms are also well studied in different 
context (8l-[l0|. Shannon entropy has also been inter- 
preted as a measure of electron correlation which is di- 
rectly related to the strength of the interaction between 
the electrons In our present study we consider in- 

teracting trapped bosons in external confinement at zero 
temperature. The gas is considered to be extremely di- 
lute as the inter-atomic separation is quite large com- 
pared to the range of interatomic interaction and the fun- 
damental features of such trapped bosons are accurately 
described by a single parameter Osc, the s-wave scatter- 
ing length. The choice of our system is quite similar to 
that of earlier studies Q , however the motivation of our 
present work is quite different, as follows. First: we are 
interested in the study of dynamics of Shannon entropies 
and the EUR in a fixed confinement but with a vary- 
ing effective interaction Nasc- For weakly interacting 
large number of bosons in the external confinement, the 
ground state properties are quite accurately explained 
by the nonlinear equation, known as time-independent 
Gross-Pitaevskii (GP) equation 



H — moj^r^ 
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g\tlj{r)\' ij(r} = fii;{r) (3) 
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where g is the interaction strength parameter given by 
g ~ 4irfi_^^jejv -g ^j^g number of atoms, m is the atomic 
mass, Use is the dimer scattering length, and lu is the trap 
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frequency. In the earher works in this direction, Eq. (3) 
has been solved and the EUR is tested for various num- 
ber of atoms 0. However in the earlier study the effect 
of nonlinearity in the dynamical feature of information 
entropies is not focussed. For this one should solve the 
time-dependent GP equation and using the dynamics of 
the many-boson wave function one can directly observe 
the effect of nonlinearity in the time evolution of infor- 
mation entropy and the EUR. Although the nonlinear 
effect in the time evolution of the condensate has been 
studied earlier [T2| - [l5j . the dynamics of information en- 
tropy and the validity of the EUR with evolution in time 
is a stronger tool to manifest the effect of nonlinearity 
and it deserves special interest. Second: we study the 
dynamics of information entropies in very tight trap. In 
the present day experiment, by controlling the external 
magnetic field, one can virtually manipulate the external 
trap size. In very tight trap, the atoms will be highly 
correlated and strongly inhomogeneous. The calculation 
of Shannon entropies will be very informative to describe 
such a highly localised and strongly inhomogeneous sys- 
tem. The study of dynamical features of information en- 
tropies and the dynamics of EUR will give rich physics. 

This paper is organised as follows. Sec. II deals with 
the brief review of time dependent GP equation and cal- 
culation of dynamics of Shannon entropies. Sec. Ill deals 
with our results. Here we analyze the dynamics of EUR 
in different trap sizes for various interaction strength pa- 
rameters g. We also present our calculated values of 
Landsberg's order parameter VL for various interaction 
strength parameters g at time t — Q. Sec. IV concludes 
with summary. 



II. DYNAMICS OF THE SHANNON 
INFORMATION ENTROPY 

We start with the zero temperature condensate in the 
harmonic trap described by the time dependent Gross- 
Pitaevskii equation as 



ih 



dt 



2m 



V 



1 



2 2 

muj r 



(4) 

To solve Eq. (4) we start with the analytic and nor- 
malized ground state solution in the absence of the non- 
linear term and propagate the GP equation with time. 
The numerical integration of the time-dependent GP is 
obtained by using a finite-difference Crank-Nicolson algo- 
rithm with a split operator technique (l6l - [l8| . We have 
verified that the numerical solution of GP equation is 
consistent with various sets of radial grids and time steps. 
In the time-dependent solution although the probability 
is time independent, i.e stationary in time, the conden- 
sate wave function has time dependent phase factor as 
exp(— i/xt), where /j, is the energy in oscillator unit. Thus 
the effect of perturbation comes through the nonlinear 
term which involves the time dependent wave function 
in the previous time step and allows us to monitor the 



condensate motion with time. 

Next we calculate the time-dependent momentum space 
wave function (l){k, t) by taking Fourier transformation of 
the time-dependent position space wave function Tp{r,t) 
at each time step as 



(27r)3/2 



-ikr 



cPf. 



(5) 



Here 'ip{r,t) is the normalized ground state solution of 
Eq. (4) and is given by ipif, t) = --i=i/)(r, i), where ip{r, t) 
is the radial part of the wave function. The angular part 
of the ground state solution is a constant and is taken 
care by the constant pre-factor along with the normaliza- 
tion. Now due to the presence of spherically symmetric 
trap, all directions are equivalent and we set the polar 
axis of the spherical polar coordinate along the direction 
of k. Then performing the integrals over the angular co- 
ordinates first, the 3Z) integral of Eq. (5) simplifies to an 
integration over the radial coordinate only. 



ilj{f,t) e 



k 



rip{r,t)sin{kr)dr. (6) 



Here k is the magnitude of k. We perform this integral 
over r by using a 32-point gaussian quadrature. We trun- 
cate the sum at r = rmax which is detrmined from the 
asymptotic behavior and normalization of ip{r,t). 

For a three-dimensional system the information en- 
tropy in the position space is calculated from the density 
distribution p{r,t) = t)p as 



Sr{t)^- / p{r,t) ln(p(r,t)) d^r 



(7) 



The corresponding information entropy in the momen- 
tum space is determined from the momentum distribu- 
tion ri,(fc, t)= as 



n{k,t) In {n{k,t)) d^k 



(8) 



If p{f,t) and n{k,t) arc normalized to unity, the joint 
entropy at time t is defined as S{t) = Sr{t) + Sk{t). At 
t = as there is no effect of nonlinear terms, S{t = 0) 
should obey the existing entropic uncertainity relation 
(EUR) [H given as 



S{t^O) =Sr{t = 0) + Sk{t^O) 
>^™n = 3(l + ln(^))c: 



6.434 



(9) 



It is already pointed out in the study of different quantum 
systems that the EUR is stronger than the Heisenberg's 
uncertainity relation due to the following reasons [6|. 
First: Heisenberg's relation can be derived from the EUR 
but its reverse is not true. Second: Heisenberg's relation 
depends on the state of the system but the EUR does 
not il. 
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Although Eq. (5) - (9) have been used earher in the 
study of interacting bosonic systems, however it needs 
additional discussion regarding the validity of the ex- 
pression for a condensate having finite size. Although 
Eq. (7) and Eq. (8) are apparently valid for single par- 
ticle only, however its extension and use for the con- 
densate is equally valid. Above the critical temperature 
Tc although there is a contribution coming from non- 
condensed atoms, but at the critical temperature Tc 
nk) there is a macroscopic occupation of the bosons. At 
such low temperature, the de Broglie wave-lengths of the 
neighbouring atoms overlap, the individual atoms lose 
their quantum identity. Thus the whole bosonic cloud 
is treated as a single quantum stuff which is well de- 
scribed by the zero-temperature GP equation. This is 
also in perfect agreement with the framework of general- 
ized Bogoliubov prescription, where the condensate wave 
function is the expectation value of the field operator. 
Thus the fluctuation in field operator is neglected. This 
is quite justified as the depletion of the condensate at 
the the temperature ~ nk is negligible. Thus Eq. (7) and 
Eq. (8) basically determines the information entropy as- 
sociated with the whole condensate and not of individual 
bosons. 



III. RESULTS 

We consider the condensate of ^^Rb atoms with asc 
= 100 ao and initially choose the trap frequency lo as 
77.87 Hz which mimics the JILA experiment [l^. With 
these parameters we solve numerically Eq. (4) for various 
number of bosons in the trap. We calculate the conden- 
sate wave function in the coordinate space with varying 
nonlinearity g. In this work we express all the quanti- 
ties in oscillator unit (o.u.). In oscillator units (o.u.) all 

lengths are expressed in units of aho (= •\/^') and all 
energies in units oifvjj, uj being the trap frequency. Since 
the oscillator units (o.u.) are frequency dependent, for 
our present work we define them {i.e. aho and hio) with 
respect to the frequency wq (=77.87 Hz) of the JILA 
trap [l^. First to check the consistency of our results 
with earlier published results Il,Sli2|, we solve the time- 
independent GP equation (Eq. (3)) with varying number 
of bosons N (or equivalently with varying strength of 
interatomic interaction g) for pure harmonic trap. In 
Fig. 1(a) we plot the condensate wave functions |-0| for 
various number of atoms in the trap. For a fixed exter- 
nal trap the wave function \'ip\ gradually expands with 
increase in number of bosons as the strength of repulsive 



0.4 

0.35 



interaction g 



) increases with increase in N . 



The momentum space wave function |0| is calculated 
by Fourier transformation of the coordinate space wave 
function \'ip\. In Fig. 1(b) we plot condensate wave func- 
tions in momentum space \4>\ for the same values of N as 
in Fig. 1(a). We see the expected reciprocal behavior be- 
tween the position and momentum space wave functions 
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FIG. 1: (color online) (a) Plot of the condensate wave func- 
tion for repulsive BEG of *^Rb atoms in harmonic trap for 
various A*', (b) Plot of the corresponding momentum space 
wave functions. 



in accordance with the Heisenberg's uncertainty princi- 
ple. Initially for small N , the effect of nonlinear inter- 
action is quite small and the condensate wave function 
is very close to gaussian. The momentum space con- 
densate wave function is also close to gaussian. As we 
increase TV, the net nonlinear interaction increases, the 
coordinate space condensate wave function |-0| spreads 
out. The corresponding momentum space wave function 
101 squeezes accordingly. For large N (=10,000) jV'l is 
quite diffuse in coordinate space and is sharply peaked 
in momentum space and has a small kink which is due 
to the presence of finite size trap. Similar kink is also 
observed in earlier studies [2^. This also nicely demon- 
strates how the presence of the harmonic trap modifies 
the momentum distribution. Next we calculate the infor- 
mation entropy Sr b.s Sr = — J p{r) In {p{r)) d^r and the 
conjugate information entropy in momentum space Sk as 
Sk = -J n{k) In {n{k)) d^k. In Fig. 2 we plot Sr and Sk 
as a function of non linear interaction g. With increase in 
particle number N, the strength of the net repulsive in- 
teraction g increases, the coordinate space wave function 
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FIG. 2: (color online) a) Plot of Sr (blue dotted curve) and 
Sk (red smooth curve) for various nonlinearity g; b) Enlarged 
view of the -ve g part. 



[■(/'(r)! delocalises and consequently the position space en- 
tropy Sr gradually increases. It signifies that the associ- 
ated uncertainty in the coordinate space increases. The 
corresponding momentum space wave function |'/>(fc)| is 
localised and the momentum space entropy Sk gradually 
decreases with increase in particle number N (i.e. g). It 
indicates that the associated uncertainty in momentum 
space decreases when net effective repulsive interaction 
increases. This inverse behavior between the position 
and momentum space entropies in Fig. 2 is interpreted 
as a consequence of the EUR and it also reflects the re- 
ciprocal behavior of wave functions in position and mo- 
mentum spaces which is the consequence of Heisenberg 
uncertainty relation. At g = 0, we have the noninteract- 
ing inhomogeneous ideal Bose gas. The important conse- 
quence of the inhomogeneity of the system is that BEC 
appears not only in momentum space but also in coordi- 
nate space. In the absence of nonlinear interatomic inter- 
action, the condensate wave function becomes perfectly 
gaussian both in coordinate and momentum spaces. Con- 
sequently the Shannon information entropy in position 
space Sr and that in momentum space Sk become equal 
to each other and we find Sr = Sk = 3.217. Thus the 



total information entropy S = Sr + Sk becomes exactly 
equal to 6.434 which is the lower bound of EUR. It is 
to be noted that the gaussian wave functions correspond 
to the coherent state with the associated Heisenberg's 
uncertainty product being minimum. For attractive po- 
tentials {g < 0), the behavior is opposite to that observed 
for repulsive potential. It is well known that in case of at- 
tractive BEC, the condensate tends to increase its density 
in the centre of the trap in order to lower its interaction 
energy. This tendency is balanced by the zero-point ki- 
netic energy which can stabilize the system. However for 
large number of bosons in the trap, the central density 
grows too much and the kinetic energy can not balance 
it anymore. The system thus collapses when the number 
of particles in the condensate exceeds a critical number 
Ncr- In the mean-field CP equation the critical number 
is determined in the following way. For attractive inter- 
action and N < Ncr, the condensate is always metastable 
and the energy functional has a local minimum. When 
N increases, the depth of local minimum decraeses and 
exactly at TV = Ncr, the minimum vanishes, and the CP 
equation has no solution. The calculated value of stabil- 
ity factor for a spherical trap is ^"^^""^ = 0.575 [l^l- In 
our present study we observe that with increase in effec- 
tive negative interaction Sr sharply falls which indicates 
that the condensate wave function in coordinate space 
\tp\ squeezes in the centre of the trap. Whereas the cor- 
responding momentum wave function |0| spreads out and 
Sk sharply increases as expected. Finally at 17 = —7.23, 
both Sr and Sk diverge. This corresponds to the stabil- 
ity factor kcr ~ ^"^^"""^ = 0.575 which is same as that 
obtained from CP t&eory. Thus we can correlate the 
collapse of the attractive BEC with the simultaneous di- 
vergence of Sr and Sk ■ The divergence in Sr nicely shows 
how the probability distribution in the coordinate space 
becomes sharply localised at the centre of the trap, like a 
(5-function. Thus our present study gives us more physi- 
cal insight in the description of the collapse for attractive 
forces. 




In N 

FIG. 3: (color online) Plot of the total information entropy 
S vs InN. Dots present the numerically calculated data and 
the straight line is the plot of Eq. (2) with calculated values 
of a and b. 
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It is an established fact that for repulsive bosons in the 
trap there is no limitation about the maximum number of 
atoms as the condensate is always stable for any number 
of atoms. Thus if we continue to increase N, Sr will 
increase and Sk will decrease accordingly as described 
earlier. However due to the inhomogeneous nature of 
the system the increase in Sr is not completely balanced 
by the corresponding decrease in Sk- Consequently the 
total entropy S will increase as shown in the Fig. 3 where 
we plot total entropy S' as a function of InA^. We find 
the plot of S vs \nN to be a straight line of the form 
S = a + blnN and it conforms to the earlier observation 
by Massen et. al. [1]. Our calculated values of a = 6.059 
and b = 0.065 are also in good agreement with earlier 
findings [1]. 

Next, we are interested in study of the dynamical 
features of entropies in both the harmonic and tight 
trap. We first calculate the time dependenct local den- 
sity in coordinate space p{f, t) and in momentum space 
t) from the corresponding condensate wave functions 
|?/'(r, t)| and respectively. Then from the local 

density p{r,t) = we calculate Sr{t) by utilizing 

Eq. (6) and from n{k,t) = we determine Sk{t) 

by Eq. (7) and then the total entropy S{t). We determine 
Sr{t), Sk{t) and S{t) with increasing interaction strength 
g to observe the effect of nonlinearity in the time evolu- 
tion of the information entropies and the dynamics of 
EUR. 

In pure harmonic trap the interacting bosons are inho- 
mogeneous. Thus the increase in Sr with particle number 
does not cancel exactly with the decrease in Sk or vice 
versa. However the effect of inhomogeneity is quite small 
for very small number of atoms in the trap. At the same 
time due to the presence of nonlinear interaction the dy- 
namical features of Sr and Sk may exhibit some nonlinear 
effect in time evolution. As the nonlinear effect increases 
with particle number, we expect to get strong nonlinear 
features even in pure harmonic trap when the effect of 
nonlinear interaction g is quite large. Thus depending 
on the total number of atoms in the external trap, the 
time evolution of Sr, Sk and the EUR can highlight the 
point of nonlinearity as the Shannon entropies are the 
best measures of localisation or delocalisation of the dis- 
tribution function. The results for g = 10 are presented 
in Fig. 4. For g=10, the number of bosons is TV = 184. 
As the system is very dilute, for such weakly interacting 
system the corresponding distribution function is close 
to gaussian.Thus we predict linear, periodic, behavior in 
Sr{t), Sk{t) and S{t). Total entropy S{t) starts from 
6.434 at t = and then exhibits regular oscillation with 
very small amplitude. The increase in Sr{t) and the cor- 
responding decrease in Sk (t) with time perfectly satisfies 
the physical meaning of the inequality Eq. (8) [l|. We 
also note that the maxima in Sr{t) corresponds to min- 
ima in Sk{t) and vice versa. Thus the joint measures of 
the uncertainty clearly signifies that for few hundreds of 
bosons, the system is very close to linear. The small effect 
of nonlinearity is smeared off by the external trap as the 



interaction energy is neglible compared to trap energy at 
zero temperature. 
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FIG. 4: (color online) Plot of Sr{t) (red smooth curve), Sk{t) 
(green dashed curve), and S{t) (blue dotted curve) with time 
t for g — 10 and uj = I.Ocjq . 

In atomic physics the total entropy maintains some rig- 
orous inequalities which can be derived using the EUR. 
Sr and Sk are fundamentally related with the total ki- 
netic energy T and the mean square radius (r^) of the 
system Q . It can be shown [2l| that 

Srr»n < S{t) < Srnax{t) , (10) 

where the lower limit of total entropy is the previously 
introduced constant and the upper limit of the total en- 
tropy is 

^™a.W = 3(l + ln(^)) + ^ln(^(r2),r(t)) , (11) 

where (r^)t is the mean square radius at time t and T{t) 
is the kinetic energy of the system. In this paper we an- 
alyze the time evolution of the inequality by calculating 
time evolution of (r'^)t and T(t) and Smax(t) according 
to Eq.(lO). We study how the above inequality Eq. (9) 
is maintained with time. 

We plot Smax{t) with total S{t) in Fig. 5. for g = 10. 
The inequality is maintained with time nicely. So the 
system is very close to linear. 

Next we observe the effect of the change of trap size on 
the dynamics of Sr, Sk and S for few hundreds bosons 
[g = 10). The trap size can be changed by varying the 
trap frequency ui (since the trap width aho)- Here we 
express various trap sizes by their corresponding trap fre- 
quncy u), measured as a multiple of the frequency of 
the JILA trap [l^. The results with uj ~ 0.85a;o (wide 
trap) and lo = 2.Quiq (very tight trap) are presented 
in Fig. 6 and Fig. 7 respectively. For uj = 0.85cjoj the 
external trap spreads compared to pure harmonic trap. 
The corresponding condensate wave function in coordi- 
nate space spreads accordingly and the momentum space 
wave function is squeezed. It is reflected in Fig. 6 nicely. 
As |'(/'(r, i)| is diffuse in such a wide trap Sr{t) increases 
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FIG. 5: (color online) Plot of S{t) (red smooth curve) and 
{t) (green dashed curve) with time t iov g = 10. 
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FIG. 6: (color online) Plot of Sr{i) (red smooth curve), Sk{t) 
(green dashed curve), and S{t) (blue dotted curve) with time 
t for g = IQ in wide trap (u) = O.SScjq)- 



and Sk{t) decreases. For such small number of bosons in 
wide trap the system is less inhomogeneous and very close 
to linear. Thus increase in Sr{t) and the corresponding 
decrease in Sk{t) almost cancell each other and the total 
entropy makes a regular periodic oscillation with small 
amplitude. However the effect of larger trap size is clearly 
seen in evolution of total entropy S{t). For uj = I.Owq 
(pure harmonic trap), we observe five distinct maxima 
in the total time period (t = to t = 8). Whereas for 
w = 0.85wo, S{t) exhibits four distinct maxima in the 
same period. The position of the first maximum in S{t) 
occurs aX t = 0.8 for lj ~ l.Owo wheras the same occurs 
in S{t) at t = 0.96 for w = 0.85cjo- Thus the wavelegth 
increases in larger trap size which is in perfect agree- 
ment with earlier observation that effect of correlation 
decreases in larger trap p^ . 

The results for tight trap with tu = 2.0cjo are presented 
in Fig. 7. In such a tight trap the system is strongly in- 
homogeneous but not strongly nonlinear as the number 
of atoms are just few hundreds. As now the condensate 
wave function is strongly localised in coordinate space, 
the corresponding uncertainty decreases and the 5*^ also 



FIG. 7: (color online) Plot of Sr{t) (red smooth curve), Sk{t) 
(green dashed curve), and S{t) (blue dotted curve) with time 
t for g — 10 in very tight trap {co — 2.0u}o). 



decreases sharply. Still as the system is close to linear 
due to such small number of atoms, Sr(t) maintains reg- 
ular periodic oscillation. Whereas we observe reciprocal 
behavior in Sk{t) which increases and exhibits periodic 
oscillation, the total entropy S{t) also shows regular pe- 
riodic behavior but with many oscillations compared to 
Ld = O.SSwo and to = I.Owq cases. The first maximum 
now appears at i = 0.32 which shows that wavelength of 
the oscillation becomes small and indicates high spatial 
coherence in tight trap [2^ . 

To observe the effect of nonlinearity we do the same 
for large number of particles which corresponds to large 
nonlinear interaction. Results for g = 100 are presented 
in Fig. 8 where we plot the evolution of total entropy 
S(t) for UJ = 0.85a;o, uj ~ I.Owq and uj = 2.0^0- For pure 
harmonic trap with N = 1840 {g = 100), we observe 
that S{t) starts with 6.434 at t = and then oscillates 
with large amplitude. Comparing with a similar figure 
for N = 184, the amplitude of oscillation is not abso- 
lutely constant throughout the whole period. It shows 
that though the system is not nonlinear, however it goes 
away from equilibrium. This is expected as the increase 
in Sr{t) does not cancel the decrease in Sk{t) as in ear- 
lier case. Thus for g = 100, the effect of inhomogcnity 
is visible in the dynamics of S(t). For w = 0.85woj we 
also observe large amplitude oscillation and wavelength 
of oscillation increases as the first maximum occurs at 
t ~ 0.64 for UJ — 1.0 whereas the same occurs at t = 0.8 
for UJ = 0.85wo. However for uj ~ 2.0a;o we observe the 
effect of inhomogcnity and nonlinearity in the dynamics 
of S(t). In Fig. 8 we observe strong nonlinear and ape- 
riodic fluctuating oscillation in S(t) which indicates that 
the system is far away from linearity. 

The results for g — 500 in various trap sizes are pre- 
sented in Fig. 9 where we observe very sharp change in 
the amplitude of oscillation in S(t) for different trap sizes. 
For very tight trap as bosons become strongly inhomoge- 
neous and highly correlated the effect of nonlinear inter- 
action becomes strong and S{t) exhibits strong aperiodic 
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(b) UJ = 1.0 




FIG. 8: (color online) Plot of S{t) with time t for 5 = 100 
and various trap widths. 



oscillation in S{t). 

In ordinary thermodynamics, entropy is a measure of 
disorder in the system. However, defining disorder A as 
A = g ^ , where S and Smax refer to the same equation 
of state, Landsberg defined a order parameter Q [2^ by 

= 1 — S I Smax- = 1 corresponds to perfect order and 
f2 = corresponds to randomness. It has been shown 
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FIG. 9: (color online) Plot of S{t) with time t for g = 500 
and vaxious trap widths. 

by a whole class of examples that the entropy S and the 
order parameter fl are two decoupled concepts 0, [2^ . 
For the sake of completeness, we have also calculated the 
Landsberg's order parameter i7 at i = 0. In Table 1 we 
present the values of f2 for various nonlinearities. 



8 



TABLE I: Landsberg's order parameter for various nonlinear- 
ity. 
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Q. 


10 


9.704 X 10~^ 


100 


9.851 X lO"'"^ 


500 


1.961 X 10""* 



It implies that by increasing the number of particles, 
the system becomes more ordered. This observation 
agrees well with earlier observation in atoms and clus- 
ters 0, Q . This is also consistent with the fact that the 
entropy and order are decoupled unlike the case in ther- 
modynamics [23j . 

Lastly our calculated values can be correlated with 
the experimental observation. As stated earlier that the 
most relevant feature of trapped Bose gas is that they 
are inhomogeneous due to external trapping potential. 
The important consequense is BEC shows up not only in 
momentum space like noninteracting homogeneous ideal 
Bose gas, but also in coordinate space [24[. This dual pos- 
sibility of investigating the effects of trapped condensa- 
tion and associated measurement of both the density and 
velocity distribution are very interesting and novel. The 
density distribution of the atoms in the trap as a func- 
tion of coordinate is measured by the optical method [l^ . 
It also provides the measurement of size of the cloud. 
Whereas other kind of distribution is the velocity distri- 
bution in momentum space which is measured by timc- 
of-flight measurements. The atoms were left to expand 
by switching off the confining trap and the velocity distri- 
bution is imaged. Thus the dual possibility to obtain the 
density distribution corresponds to the Fourier transfor- 
mation of the condensate wave function in the momen- 
tum space. As the density in coordinate and momentum 
space are the basic ingredients to calculate the Shannon 
information entropy, our theoretical results can be cor- 
related with the experimental observation. This is also 
nicely demonstrated in Fig. 2(a) and 2(b), where the 
opposite diverging behavior in Sr and Sk nicely demon- 
strate the experimentally observed critical point of sta- 
bility. Again the maximum value of information entropy 
can be accurately determined from Eq. (11) where the 
size of the condensate < >t and the kinetic energy 
of the system are known from the experimental results. 
Thus our theoretical results provide rich physics which 
can be tested in modern experiments. 

Another relevant question is how to correlate the the- 
oretical T = temperature results with those of ex- 
perimentally measurable quantities which would be for 
T > 0. In the real experimental situation the condensate 
temperature is ~ nk as said earlier. Thus in the pres- 
ence of thermal cloud, there is an interaction between 
the condensate and non-condensed atoms and thermal 
fluctuation also plays an important role. However in the 
present study we consider the T>0 but T<Tc (Tc being 
the critical temperature) where there is a macroscopic 



population in the single lowest quantum state, i.e. the 
condensate is absolutely formed and there is no thermal 
fluctuation. Thus our system does not have thermally 
excited atoms. Apparently it may contradict the exper- 
imental situation, however in the presence of external 
trap as the interaction energy is very small compared to 
the trap energy, one can safely use the zero-temperature 
GP equation in the description of condensate. However 
when T>Tc, there will be thermal distribution of atoms 
in the higher states, the ground state is depleted and to 
study such system we must use temperature dependent 
formalism. Thus our T = temperature results are par- 
ticularly valid in the length scales from the interparticle 
distance to the healing length which is the relevant range 
in current experiments where na^c ^ 1- 



IV. CONCLUSION 

Shannon entropies are generally used to examine the 
localisation of the distribution function in the position 
and momentum space. In the present paper the dynam- 
ics of information entropies and the EUR are studied with 
varying interparticle potential of the interacting trapped 
bosons and varying trap width. The system draws spe- 
cial attention due to the presence of external trap. We 
observed that for weak interaction (i.e small g) the sys- 
tem shows periodic behavior irrespective of trap widths. 
However with increasing interaction strength g, irregu- 
lar behavior comes in and also the effect of varying trap 
width becomes more pronounced. We stress that since 
we choose the parameters like ujq, cisc etc. as those of 
JILA experiment and also restricted N up to few thou- 
sands only, GP equation is good enough to describe it. 
For more dense system with stronger interaction strength 
(i.e. larger N and Qsc) inter-atomic correlations become 
important. So some other theoretical technique incorpo- 
rating correlations needs to be applied. Another impor- 
tant observation is that the EUR is maintained through- 
out time for all cases and thereby reaffirms its generality 
as an earlier study has verified it numerically for differ- 
ent number of bosons in the trap. The increase of f2 with 
g (i.e. N, as we kept age fixed) is consistent with ear- 
lier study and confirms again that the entropy and order 
are decoupled. However the dynamic evolution of fl with 
time for various g remains to be studied. Finally as de- 
scribed earlier, all our results are directly related to the 
experimentally measurable quantities and therefore they 
are experimentally accessible even with available set up. 
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